Short Range Scaling Laws of Quantum Gases With Contact Interactions 
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The probability amplitude for A'' particles in a quantum gas with negligible range of interparticle 
interaction potentials to come to a small region of size r scales like r^. It is shown that 7 is 
quantitatively related to the ground state energy of these A'^ fermions in the unitarity limit, confined 
by an isotropic harmonic potential. For large N, the short range density distribution of these 
particles is predominantly the same as the Thomas-Fermi profile of the gas in the unitarity limit 
confined by such a harmonic potential. These results may shed light on strongly interacting ultracold 
atomic Fermi gases, in a trap or on an optical lattice. 
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The recent realization of degenerate atomic Fermi 
gases with large scattering lengths, near magnetic-field 
modulated Feshbach resonances P, 0i S Hfi provides 
motivation to study a simple model of quantum gases, 
in which the interaction potentials are of zero range, and 
the scattering amplitude between two interacting parti- 
cles is characterized by a single parameter, the s-wave 
scattering length a ,6j . In reality, if the range of inter- 
action potential between particles is much shorter than 
the other length scales, the scattering length a and the 
typical de Broglie wave length of particles, such a simple 
model should be justifiable. 

An accurate theory about this model must respect all 
the features of the many-body wave functions of this sys- 
tem. So the structure of such wave functions deserves 
in-depth investigations. 

In this paper we study an aspect of the short-range 
structure of these wave functions. How does the full 
many-body wave function ip (regardless of ground state 
or any excited state which is present in the density oper- 
ator at finite temperatures) scale with r, when N of the 
particles come to a small spatial region of size r ^ 0? 
Here r is still much larger than the actual ran_ge of the 
interaction potential. We expect a scaling law 
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or, more intuitively (but less rigorously), tp ^ r'^ . The 
partial differential in Eq. is defined as follows: 1) 
reduce the distance between each of the N particles and 
a fixed reference point O by an amount proportional to 
such a distance, without changing the particle's direction 
with respect to O; 2) hold all the other particles at their 
locations, r is any length scale characterizing the spatial 
extension of these N particles. 

A major new result in this paper is that 7 has an in- 
teresting asymptotic formula at large N, and it is deeply 
related to the ground state energy of the same quantum 
gas in the unitarity limit in which \a\/l goes to infin- 
ity. / is the average distance between particles. A broader 



new finding is that for arbitrary N, 7 is quantitatively 
related to the ground state energy of these N particles 
in the unitarity limit, confined by an isotropic harmonic 
potential. 

7 is a function of and usually it does not depend 
on the relative configuration of these N particles. 7 does 
not depend on the scattering length a, because when two 
scattering particles come to a distance r ^ |a|, the wave 
function is dominated by the term proportional to 1/r. 
7 does not depend on the temperature of the system, 
either, because for any energy eigenstate, the typical de 
Broglie wave length of a particle is very long compared 
to the small value of r. 

In the case of a two-species Fermi gas, with only an in- 
terspecies s-wave contact interaction and a finite scatter- 
ing length, 7(1,1) = -1 trivially 9']. Here j{Ni,N2) = 
j{N2,Ni) stands for the value of 7 for A^i fermions of 
one species and A^2 ones of the other species coming to a 
small region of space. 

7(1,2) was calculated in f-^d its value is 

—0.2272757 • ■ • . This is associtated with a p-wave scat- 
tering of three fermions, one of species A, and the 
other two of species B 10!|. If the total orbital angu- 
lar momentum of these three fermions is precisely zero, 
7 = 0.1662 • • • jlfl|. For a general wave function ip, this 
greater value of 7 only appears at a subset of the config- 
uration space, with measure zero [l l| . From this point 
on, we will only consider the least possible 7 in any case. 

It is feasible, but nontrivial, to calculate 7(-/Vi, A^2) for 
all other values of A^i and A^2- 

Here we study the asymptotic behavior of '-f{N/2, N/2) 
at large N. We expect, before we do any concrete cal- 
culations, that this asymptotic behavior must be closely 
related to the ground state of the 2-component Fermi 
gas in the unitarity limit. The reasons are: 1) when N 
fermions come to a small region in space, a can be re- 
garded as infinity, also the total energy of these fermions 
is negligible, compared to the energy scale set by Heisen- 
berg's uncertainty principle h'^/mr'^, where m is the mass 
of each particle and h Plank's constant divided by 27r, 
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2) for large N, this "few" -body system should approach 
the many-body system, and 3) the least possible 7 cor- 
responds to the ground state, as will be clear shortly. 

For a little more generality, we consider a g-component 
Fermi gas of particles of identical mass to, with either 
contact interactions or no interactions at all, and approx- 
imately calculate the value of 7 for N fermions, with N/g 
of each species, coming to a small regkm of space. In a 
hyperspherical coordinates formalism these fermions 
can be described by a hyperradius R, and a set of hyper- 
angles fl. 

R' = j:rl (2) 

i=l 

where is the distance between the i-th fermion and 
the center-of-mass of the N fermions. The Schrodinger 
equation for these N fermions is 

where the operator K is the hyperangular contribution 
to the total energy. We have omitted the mass-reduction 
effect because of the large N, and omitted the small con- 
stant term —4 in the numerator of the fraction 3N/R. 
We also omit the total energy of the N fermions at small 
-R, so the hyperradial part in the Schrodinger equation 
cancels the hyperangular part. 

We need to estimate K = {Kip)/ip, in order to es- 
timate 7. Physically, K should be the lowest possible 
total energy of the N fermions at a given R, because 
higher energies lead to greater values of 7 which will be 
overwhelmed by the least 7 at small distances. So we 
should use the formula for the ground state energy of 
the Fermi gas in the unitarity limit. In the local den- 
sity approximaion 13], which should be valid when N is 
large, the ground state energy per particle is dK/dN = 
(3/5)(l + I3)h'^kl/2m, where kp = {Gn^p/g)'^/^ is the 
Fermi wave vector, p is the local number density of 
fermions, and (3 a universal constant 14»|, dependent only 
of the composition of the Fermi gas in the unitarity limit. 
Obviously the density profile of the N fermions should be 
isotropic, if we are to achieve the minimum energy at a 
given R. So p ~ pir), where r is the distance to the 
center-of-mass. 

We now have three equations, for the number of scat- 
tering fermions N, the hyperradius R, and the energy K 
at this given hyperradius. 

/>oo 

N = p(r)47rr^dr, (4a) 
^0 

/>oo 

i?2 ^ / r^p{r)4:nr^dr, (4b) 
Jo 

r'Jl±J^('-^Y\ir)^fH.r^dr. (4c) 
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To minimize K at given values of A'" and R^ , we use two 
Lagrange multipliers, fj, and v, and minimize K — fiN — 
vR^ by adjusting p{r). The only physically meaningful 
solution to this problem is 

for r < ro, and p{r) = for r > r^. Here ro — \J [ij (^—v) 
and < 0. The cloud of fermions at the intermediate 
stage of scattering (with small instantaneous hyperra- 
dius i?) thus have a density distribution equivalent to 
the Thomas-Fermi profile of a Fermi gas in the unitarity 
limit confined by an isotropic harmonic potential. 

Expressing p. and ro in terms of N and i?^, and subsi- 
tuting the results into the equation for we get 

K = .g"2/^(l -f fi){%Nfl'^h^ l(X2&mR^). (6) 

Note that K is inversely proportional to the hyperradius 
squared, as a consequence of the unitarity limit |l4j. As- 
suming ij] <y: R'^ we then derive from Eq. ^ that 

7(7 - 1 + iN) = .g~^/^(l -t- l3){QNf'^ l&A. 

Pauli blocking suppresses the amplitude for many iden- 
tical fermions to come to a small region of space, so only 
the positive solution should be retained, and 

7 = g-^'^'^/YTpmf'^l^ - 3N/2, (7) 

with a relative error likely to be of the order A^~2/3 
large A^. 

In the case of the 2-component Fermi gas with inter- 
species contact interaction, 1 + f] ^ 0.44 according to the 
latest numerical calculations 16], so 

7(A^/2, A^/2) < 0.72A^'^/^ - 3A^/2, for large N. (8) 

Equation Q also holds for the single-component Fermi 
gas with no interactions, for which 5 = 1 and /3 = 0. In 
this case 7(A^) can be exactly determined for any A^, as 
follows. 

For A^ = 1, the wave function is in most regions a 
nonzero value, so 7(1) = 0. For A^ = 2, we may consider 
a Slater determinant of two base functions in a small 
region of space, 1 and x, and the wave function is of the 
form xi — X2, where the subscripts label the fermions. 
So 7(2) = 1. Actually there are three independent linear 
functions of r, namely x, y, and z. We may "fill in" 
one more fermion, and form a Slater determinant of, say, 
1, X, and y. This means that 7(3) — 2. And similarly 
7(4) = 3, for which the constant "orbital" and the three 
linear "orbitals" are all filled by fermions. Note that this 
kind of short-range analysis is valid in virtually all regions 
of the many-fermion configuration space, except regions 
of measure zero. 
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FIG. 1: The short range scahng exponent 7 vs. the number 
of particles A*', for a single-component Fermi gas with no in- 
teractions. Black dots: the exact 7 values. Solid line: the 
asymptotic formula derived in this paper for large N, Eq. Q. 



For five fermions, the constant orbital and the three 
linear ones are not sufficient. The fifth fermion has to 
be filled into a quadratic orbital, say x'^. So the Slater 
determinant is at least a fifth order polynomial of the 
coordinates, and 7(6) = 5. There are six independent 
quadratic orbitals, x^, y^, z^, xy, yz, zx, and we can fill 
them one by one. Each time a more fermion is filled into 
these orbitals, 7 increases by 2. So 7(6) = 7, 7(7) = 9, 
. . . , and 7(10) = 15. 

This analysis can be easily generalized to arbitrary N. 
In general, either exactly all the orbitals of orders no 
larger than n are filled, ie, we have a short range closed 
shell structure, or a portion of the {n + l)-th order or- 
bitals are also filled, ie, we have a short range open shell 
structure. In the former case, N — (n-|-l)(7i + 2)(n-t-3)/6 
and 7(A^) — n{n -|- l)(n -|- 2)(n + 3) /8. In the latter case, 
N ^ {n + l)(n -t- 2)(n -I- 3)/6 + S, -f{N) = n{n + l)(n -|- 
2){n + S)/8+{n+l)S, and 1 (5 < {n + 2){n + 3)/2. In 
both cases, 7 — 

(67V)4/3/8 _ 3^/2 with relative errors of 
the order iV~^/'^ at large N. This is completely consis- 
tent with Eq. O. See Fig.jTlfor the comparison between 
the asymptotic formula and the exact 7 values. 

It is straightforward to verify the validity of Eq. JT)) for 
the g-component Fermi gas with no interactions. 

There is no reason to believe that the favorable com- 
parison between our asymptotic formula and the actual 
7's should break down for the two-component Fermi gas 
with interspecies long scattering length, although no one 
has ever computed the actual 7's for iV ^ 4. Here we 
just show that for N = 2, Eq. {Tj) gives an estimate of 
7 = —1.19, according to the latest l+P value 0|. This is 
already surprisingly close to the actual 7 value, —1. Note 
that for a 2-component Fermi gas with no interactions, 
our asymptotic formula gives 7(1, 1) = —0.27, while the 
actual 7(1,1) = 0. The magnitudes of errors are appar- 
ently similar for these two different systems. If we solve 



Eq. ^ for 1 -t- f3's from the known exact 7(1, 1) values, 
we get 1 + P = 0.538 for the 2-component Fermi gas with 
interspecies long scattering length, and 1 + fi = 1.21 for 
the 2-component Fermi gas with no interactions. Both 
of these two estimates are larger than the actual 1 + /3 
values by just about 21-22%, at such a small number of 
fermions as iV = 2. 

We may generalize Eq. lO to the case in which the 
numbers of fermions coming to a small spatial region are 
not equal among different species, but have some fixed 
ratio. In this case, we expect to approach the ground 
state of the Fermi gas in the unitarity limit with that 
ratio among the populations of different species, and will 
get some different p. Also we can study the case in which 
different components have different masses. 

We may also generalize the same analysis to other 
quantum gases with contact interactions, as long as those 
gases are stable. A prominant example is a Fermi-Bose 
mixture, with no interactions between the bosons. If 
the particles all have the same mass, and the number 
of colliding particles is evenly distributed among all the 
species, the formula for 7 is formally the same as Eq. Q, 
except that the value of (3 is different. 

Our results are a demonstration of quantum paral- 
lelism in the context of quantum gases with long scatter- 
ing lengths. Normally one regards the regimes in which 
kpa 0(1) as distinct from the unitarity limit, in which 
kpa ~ 00. We have shown, however, that such a simple 
picture is not complete. In fact, even for finite kpa, small 
portions of particles still have finite, albeit very small, 
probability amplitudes of clustering for short times, and 
when they do so, their behavior is very similar to the be- 
havior of the entire quantum gas in the unitarity limit. 
We may call this peculiar behavior instantaneous unitar- 
ity limit. Note that such a phenomenon coexists with 
other more familiar ones in any many-body wave func- 
tions, regardless of the ground state or excited states. 
Non of the many-body theories built to date reflect this 
feature. 

Our short-range structure analysis of the many-body 
wave function may shed light on the theory of the quan- 
tum gases away from the unitarity limits. One interest- 
ing question is: how to construct an accurate theory of 
a 2-component Fermi gas in the pair-wise Bose-Enstein 
condensation (BEC) states, when kpa evolves from small 
positive numbers to values of order unity? In a low- 
density expansion in the BEC regime, one is forced to 
take progressively higher orders of correlations into ac- 
count, with more and more dimers colliding at distances 
of the order a. One has to evaluate all these scatter- 
ing amplitudes, when the density is not extremely low. 
While these amplitudes may be easier to calculate di- 
rectly when the number of colliding dimers is not large, at 
large numbers of dimers, one will find the idea of instan- 
taneous unitarity limit, or, more generally, instantaneous 
high density, helpful. 
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Another possible application of the ideas in this pa- 
per is: we may refine our asymptotic formula, by taking 
all the correction terms to the Thomas-Fermi approxima- 
tion into account; then, by solving a quantum mechanical 
problem of a small number of fermions, we could accu- 
rately determine the physical properties of the gas in the 
unitarity limit. 

In the end of this paper, it should be pointed out 
that there is an exact short range scaling theorem, which 
places Eq. Q on a broader and more rigorous ground. 

Theorem. The ground state energy of N fermions of 
identical mass m with contact interactions of infinite 
scattering lengths or no interactions, confined in an 
isotropic harmonic trap of angular frequency to, is 

E^{j + 3N/2)nuj (exact), (9) 

where 7 is the short range scaling exponent defined pre- 
viously. 

Proof. Let us start with the wave function of N fermions 
with contact interactions of finite scattering lengths or 
no interactions, ip{ri,--- ,rjv)- When ^ but the 
distance between any two interacting fermions remains 
nonzero, ^p satisfies the scaling law and the Schrodinger 
equation 

JV 

^r, •V,V = 7V', (10a) 
1=1 

tfj also must satisfy appropriate boundary conditions 
whenever the distance r between two fermions vanishes, 
and is like c/r + 0{r) for a/r — > 00, or c+0(r) for a ~ 0, 
where c is independent of r. We now extrapolate the iV- 
fermion short-distance wave function tp to finite distances 
according to the scaling law, and then define 

^'^^exp(^-^f]rf^, (11) 

which clearly has the same symmetry and boundary 
conditions as ip. It is straightforward to derive from 
Eq. (|10|) that for nonvanishing distances between inter- 
acting fermions 

E + 2™""'^'' ) ^' = (7 + ^N/2)hw^P'. (12) 

So for the lowest possible 7, ■0' is the ground state wave 
function of the N fermions (with zero or infinite scatter- 
ing lengths) in the harmonic trap, with the energy given 
by Eq. ®. □ 



At large N, we can rederive Eq. by using Eq. ® 
and the Thomas- Fermi approximation ^3 of ^ ^3 • We 
see that the term —3N/2 in Eq. Q is deeply linked to 
the zero-point energies of harmonic oscillators. 
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